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Abstract 

This article presents a power-spectrum analysis of 2,350 measurements of the ^"Sr/^^Y decay process acquired over the interval 4 
August 2002 to 6 February 2009 at the Lomonosov Moscow State University (LMSU). As we have found for other long sequences 
of decay measurements, the power spectrum is dominated by a very strong annual oscillation. However, we also find a set of low- 
frequency peaks, ranging from 0.26 year"' to 3.98 year"', which are very similar to an array of peaks in a power spectrum formed 
from Mt Wilson solar diameter measurements. The Mt Wilson measurements have been interpreted in terms of r-mode oscillations 
in a region where the sidereal rotation frequency is 12.08 year"'. We find that the LMSU measurements may also be attributed to 
the same type of r-mode oscillations in a solar region with the same sidereal rotation frequency. We propose that these oscillations 
occur in an inner tachocline that separates the radiative zone from a more slowly rotating solar core. 
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1. Introduction 

Analyses of decay rates of radioactive elements acquired at 
Time Nature Explorations at the Lomonosov Moscow State 
University (LMSU) have revealed evidence for variability, 
specifically oscillations with periods of one year and of approx- 
imately one month |1]. We have noted |2] that an annual peri- 
odicity is also exhibited by data acquired at the Brookhaven 
National Laboratory (BNL) |3] and by data acquired at the 
Physikalisch-Technische Bundesanstalt (PTB) \^. Power spec- 
trum analysis of the BNL and PTB datasets also reveals a mod- 
ulation with a period of order one month, which may be due to 
solar rotation ||5H7|]. It is also probable that the modifications of 
decay rates of radioactive nuclides are influenced by the relic 
neutrino flux |8]. 

' This article presents an analysis of an approximately daily 
compilation of 2,350 measurements of the ^'^Sr/^'^Y decay pro- 
cess acquired at the Institute for Time Nature Explorations at 
the Lomonosov Moscow State University over the time interval 
August 4 2002 (2002.59) to February 6 2009 (2009.10) ^. In 
order to obtain a smoothly running set of dates that are close 
to calendar dates, we count days from 1970 January 1 as day 1 
and then convert to years by dividing by 365.2564 and adding 
1970. (This procedure has proved convenient for the analysis 
of neutrino data, which all fall into this timeframe.) 
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The possibility that the solar core may rotate more slowly 
than the radiative zone (that has a sidereal rotation rate of about 
13.5 year"') has led us to conjecture that there may be an "inner 
tachocline," separating the core from the radiative zone lliol] . 
By analogy with the "outer tachocline" that separates the radia- 
tive zone from the convection zone, this would be a localized 
region where there is a sharp radial gradient in the rotation rate. 

The well known Rieger oscillation [1 1] may be interpreted as 
an r-mode oscillation [12] in the outer tachocline [Il3|. These 
oscillations are retrograde waves in a rotating fluid with fre- 
quencies determined by the sidereal rotation frequency, vr, and 
the values of / and m, two of the three spherical harmonic in- 
dices. The values of I and m are restricted by / = 2, 3, . . . , m = 
1 , 2, . . . , Z. The frequencies relevant to processes within the Sun 
are given by 



v(l,m) — 



2mvR 
Kl+l) 



(1) 



The Rieger oscillation, with frequency 2.37 year"' correspond- 
ing to a period of 154 days, may be interpreted as an r-mode 
oscillation with / = 3, ;« = 1, where the sidereal rotation fre- 
quency is 14.2 year"'. This rotation frequency identifies the 
source of these oscillations as the outer tachocline [14]. 

We have recently carried out a power-spectrum analysis of a 
sequence of 39,024 measurements of the solar diameter made 
at the Mount Wilson Solar Observatory [15]. This analysis 
yields strong evidence for an array of oscillations that may 
be attributed to a set of r-mode frequencies (for m = 1 , / = 
2, 3, 4, 5, 6, 7, 8, 10) originating in an internal location where the 
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sidereal rotation frequency is 12.08 year '. This estimate also 
is suggestive of a source in the inner tachocline. 

The present article is concerned with measurements of the 
^"Sr/^^'Y decay process acquired by Time Nature Explorations 
at the Lomonosov Moscow State University (LMSU). Parkho- 
mov has found that these measurements reveal evidence for two 
oscillations, one with a period of one year and the other with a 
period of approximately one month, and suggests that the an- 
nual modulation may be due in part to a cosmic-neutrino com- 
ponent of dark matter (jsl 191] . 

We carry out a power-spectrum analysis of the LMSU dataset 
in Section 121 We search the power spectrum for evidence of r- 
mode oscillations in Section [3j and we discuss the results in 
Section|4] 

2. Power Spectrum Analysis 

We have determined, and corrected for, the mean decay rate 
of the LMSU measurements of the ^"Sr/^^'Y decay process, 
and then normalized the corrected measurements to mean value 
unity. The resulting normalized measurements are shown as a 
function of time in Fig. [T] The standard deviation of these data 
is 0.11%. 
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Figure 1: "'Sr/"'Y decay measurements, after adjusting for mean decay rate 
and normalizing. 

We next carried out a power-spectrum analysis of these nor- 
malized measurements by a likelihood procedure [Uct that is 
an extension of the Lomb-Scargle procedure |jl7l |l8[. The re- 
sulting power spectrum, restricted to the frequency range 0-5 
year"', is shown in Fig. |2l The power spectrum is dominated 
by a very strong peak at 1 year"', as is typical of decay-rate 
measurements [5-7]. 

If the dataset comprised both measurements and error esti- 
mates, we could examine the robustness of the power spectrum 
by sampling the data many times, drawing each measurement 
from a normal distribution with the appropriate mean and stan- 
dard deviation. We do not have error estimates for the dataset 
we are examining, but we may examine the robustness of the 
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Figure 2: Power spectrum analysis of tlie normalized data. The peak at 1 year"' 
has power 65 1 . 



power spectrum by using the bootstrap procedure | ll9j . Given 
A^ time-measurement pairs, we may extract one pair randomly, 
then draw another pair randomly, having replaced the first pair, 
and so on, until one has constructed a simulation of the A^ time- 
measurement pairs. We then form the power spectrum of this 
simulated dataset. We have carried out 100 such bootstrap sim- 
ulations, from which we form the mean power and the standard 
deviation of the power at each frequency. The result is shown in 
Fig. |3] from which it is clear that the uncertainty in the power 
spectrum, estimated in this way, is only a minor factor in our 
analysis. We therefore ignore this uncertainty from now on. 



700 




- 


600 




- 


500 




- 




1 400 

CL 




- 


300 




- 


200 




- 


100 
n 


mJ 


k^A ^JX^.^^ 



12 3 4 5 

Frequency (year"^) 

Figure 3: Power spectrum formed from 100 bootstrap simulations of the nor- 
malized data. Mean value in red, mean ± standard deviation in blue. 



3. Evidence for r-Mode Oscillations 

We now seek to determine whether there is any value of 
the rotation frequency v« for which the power spectrum shows 
peaks corresponding to values of the r-mode frequencies given 
by Eq. [1] This requires the simultaneous examination of the 
power at several different frequencies. A convenient procedure 
is to use the "Combined Power Statistic" 12011 . defined as fol- 
lows: If U is the sum of n powers, 



t/ = 5i +S2 + ■.■ + S„, 
then the CPS is defined by 

G(t/) = f/-ln(l + f/-t 



:U' 



1 



in - 1)! 



:U 



n-\ 



(2) 



(3) 



This statistic has the following property: If each power is drawn 
from an exponential distribution, as is appropriate if the power 
is formed from a time series dominated by normally distributed 
random noise [18] 18], then G also conforms to an exponential 
distribution. 

We show this statistic in Fig |4] for vr in the range 10- 
15 year"' (which is the search band we have adopted, as 
a convention, for internal solar rotation rates), combining 
power measurements for the following r-modes: m - \,l - 
2, 3, 4, 5, 6, 7, 8, 9. We see that there is a well-defined peak close 
to vs=12.0year"'. 
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Figure 4: The Combined Power Statistic, formed form m = 1,1 = 
2,3,4,5,6,7,8,9, r-mode frequencies. Tlie peak is found at v=12.08 year"', 
with CPS value 215.6. 

Figure|5]shows the power spectrum over the frequency range 
0-5 year"', with arrows indicating the r-mode frequencies for 
m = 1,/ = 2,3,4,5,7,9, and for vr=12.08 year"', which is the 
sidereal rotation rate derived from our analysis of Mt Wilson di- 
ameter data 115]. We see that there is good agreement between 
these values and peaks in the power spectrum. (There are no 
peaks corresponding to m = 1, / = 6 and 8.) 
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Figure 5: The Combined Power Statistic, formed form m = 1,1 = 
2,3,4,5,6,7,8,9, r-mode frequencies. The peak is found at y=12.08 year"', 
with CPS value 219.5. 



As a test of the significance of this result, we have carried 
out 10,000 shuffle simulations ]21] of the data, with the result 
shown in histogram form in Fig. |6] None is larger than 12.8, 
whereas the actual value is 215.6. 
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Figure 6: Histogram of the combined power statistic fornied from 10,000 
shuffle-simulations for v=12.08 year"'. None is lai'ger than 12.8, whereas the 
actual value is 215.6. 



4. Discussion 

We show in Fig. Q the power spectrum recently derived 
from Mt Wilson diameter measurements lll5ll . This is remark- 
ably similar to the power spectrum derived from ^"Sr/^*'Y data, 
shown in Fig. |5] suggesting that they have a common origin. 
Each of these datasets exhibits periodicities that may be at- 
tributed to a complex of r-mode oscillations. The Mt Wilson 



data show evidence of oscillations corresponding to m = 1,1 - 
2,3,4,5,6,7,8, 10 and a sidereal rotation frequency of 12.08 
year"', the peak at / = 9 being missing. By contrast, we see 
from Fig. |5]that, in the LMSU power spectrum, the peaks at 
/ - 6 and 8 are missing. These missing modes may not be 
significant, since studies of the Rieger [11] and similar oscil- 
lations show that solar r-modes oscillations are quite typically 
intermittent. 
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Figure 7: Power spectrum formed from Mt Wilson diameter data, showing the 
locations of peaks expected for r-modes for sidereal rotation frequency 12.08 
year"'. Red an'ows indicate m = 1,/ = 2, 3,4, 5,7, 8, 10. 

As mentioned in the Introduction, the Rieger and related os- 
cillations may be attributed to r-mode oscillations that develop 
in a region where the sidereal rotation rate is near 13.7 year"', 
placing it in the familiar tachocline 1141 . This suggests that r- 
mode oscillations are unstable and can grow to large amplitudes 
in regions where there is a radial gradient in angular velocity. 
(In this connection, it is interesting to note that Zaqarashvili et 
al. have shown theoretically that a latitudinal gradient in an- 
gular velocity can lead to instability [22].) Since the r-mode 
oscillations manifested in the LMSU and Mt Wilson data seem 
to occur where the sidereal rotation rate is close to 12 year"' 
(much less than the rotation rate of the radiative zone), these 
results are compatible with the proposal that these oscillations 
develop in an inner tachocline that separates the core from the 
radiative zone. 

We have recently presented the results of power-spectrum 
analyses of nuclear decay data acquired at the Brookhaven 
National Laboratory (BNL) ^ M and at the Physikalisch- 
Technische Bundesanstalt laboratory (PTB) llaLZl. Analysis of 
BNL data yields evidence of periodicities at or near 1 1 .2 year" ' 
and 11.9 year"' [5]. When the BNL analysis is combined with 
an analysis of PTB data, we find the dominant oscillation to be 
at about 11.2 year"' |7j]. If this is interpreted as a rotational 
modulation, we see that these oscillations originate in a region 
with a sidereal rotation rate of about 12.2 year"', not far from 
the estimates of rotation rates inferred from our analysis of r- 
mode oscillations. Fig. [8] shows the power spectrum formed 



from LMSU data for the frequency range 10-15 year"' We see 
that the biggest peak in this range is one at 10.92 year"' with 
power S - 4.93. If interpreted as a rotational modulation, it 
corresponds to a sidereal rotation frequency of 11.92 year"', 
not far from the value (12 year"') inferred from our r-mode 
analysis. 
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Figure 8: Power spectrum formed from LMSU data for the frequency range 
10-15 year"'. The biggest peak is found at 10.92 year"', with power S =4.93. 

Concerning internal rotation rates, it is interesting to re- 
view the results of power spectrum analyses of the Super- 
Kamiokande solar neutrino data [|23|[ . The most prominent peak 
in a power spectrum analysis of Super-Kamiokande data (in 5- 
day bins) is found at 9.43 year"' [20[. If this peak is attributed 
to rotational modulation, it indicates that the Sun contains a 
region with a sidereal rotation rate of 10.43 year"', which is 
suggestive of a slowly rotating core. 

In discussing the power spectrum of Super-Kamiokande data, 
it is important to note that if these data are analyzed by the 
Lomb-Scargle procedure (as in the analysis carried out by the 
Super-Kamiokande Consortium 12411 '). this peak does not appear 
to be statistically significant. However, if the data are analyzed 
by a likelihood procedure that takes account of the error esti- 
mates and of the specific durations of the time-bins (which the 
Lomb-Scargle procedure does not), the peak is found to be quite 
significant OOll . 

If the inner tachocline extends from the core, with a rotation 
rate of about 10.4 year"', to the radiative zone, with a rotation 
rate of about 13.7 year"', we see that the rotation rate has the 
value 12 year"' in the center of the tachocline. There is a corre- 
spondence between this result and the fact that the Rieger and 
related oscillations seem to have their origin in the center of the 
outer tacholine. 
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